Reentrance effect in the lane formation of driven colloids 
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Recently it has been shown that a strongly interacting colloidal mixture consisting of oppositely 
driven particles, undergoes a nonequilibrium transition towards lane formation provided the driving 
strength exceeds a threshold value. We predict here a reentrance effect in lane formation: for fixed 
high driving force and increasing particle densities, there is first a transition towards lane formation 
which is followed by another transition back to a state with no lanes. Our result is obtained both by 
Brownian dynamics computer simulations and by a phenomenological dynamical density functional 
theory. 
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If a binary colloidal mixture is exposed to a constant 
external force which acts differently on the two parti- 
cle species, the system achieves a nonequilibrium steady- 
state in general For strong external forces, par- 
ticles of the same species, which are driven alike, align 
behind each other ^ ^] . This "slip-stream" effect results 
in macroscopic lanes spanning the whole system size pro- 
vided a threshold in the strength of the driving force is 
exceeded. Lane formation can be classified as a nonequi- 
librium first-order phase transition. The transition to- 
wards lane formation is general and occurs also in sys- 
tems different from colloids, such as dusty plasma par- 
ticles ions migrating within two-dimensional mem- 
branes 7] , granular matter 8, 9], and collective dynamics 
of pedestrians in pedestrian zones . 

A quantitative investigation of lane formation has been 
performed by Brownian dynamics computer simulations 
of driven two-component colloidal mixture in two and 
three dimensions P| . These simulation studies have been 
extended to asymmetric mixtures relatively tilted 
driving fields acting on the different particle species ^3 
and crystalline mixtures driven against each other |l2j |. 
Lane formation also occurs as a transient state during 
the two-dimensional Rayleigh- Taylor interfacial instabil- 
ity provided the interfacial tensions are small with re- 
spect to the driving force • Moreover a theory for the 
transition towards lane formation was recently proposed 
by us [31 ■ It is based on dynamical density functional 
theory, which takes the repulsive interparticle interaction 
explicitly into account, and is supplemented by a phe- 
nomenological version of the transversal particle current 
induced by the external drive. An instability from a ho- 
mogeneous state towards a laned state was predicted but 
the actual predictions for the threshold of the external 
drive were significantly smaller than in the simulation Q . 



The previous studies in Refs. 0,0] were restricted to 
high densities of the colloidal particles. In this brief re- 
port, we illustrate our observations on the low density 
regime as well to bring out the entire topology of the 
laning phase diagram. We show that there is a reentrant 
effect in driven colloidal mixtures. Reentrance occurs as 
a function of particle density at fixed strength of the ex- 
ternal drive. For increasing density and sufficiently high 
driving force, we find the sequence: no lanes/lancs/no 
lanes, i.e. lane formation occurs only in a finite density 
window. There is a simple intuitive reason for this reen- 
trance effect: for very small densities, the particles are al- 
most non-interacting, hence the drive will not induce any 
instability so that the system will stay mixed for entropy 
reasons. For intermediate densities, there is an interac- 
tion and the strong drive will put particles into lanes. For 
very high densities, on the other hand, the interaction is 
very strong exceeding the external drive so that the ex- 
ternal drive will again not induce lane formation. Our 
results are based on Brownian dynamics computer simu- 
lations. We also adopt the dynamical density functional 
theory of Ref. 13. However, we improve the phenomeno- 
logical expression for the transverse current in this theory 
and show that the modified theory captures the reentrant 
feature and brings about much better agreement with the 
simulation data for the location of the laning transition. 

The Brownian dynamics simulation model in this work 
is the same as used in our previous studies 4, 14]. Let 
us briefly summarize the essential features: we consider 
an equimolar binary colloidal system with 2N = 500 
particles in a square cell of length ^, having periodic 
boundary conditions with a fixed total number density 
po = 2N/f' = 2p\ = 2p° and a fixed temperature T. 
The effective pair potential between two colloidal parti- 
cles at a separation r is modeled as a screened Coulomb 
interaction 2'] 



V{r) — VqU exp [—K{r — a)/a\/ r. 



(1) 
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where Vq is an energy scale, a is the length scale, defining 
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FIG. 1: Simulation snapshots for a fixed external force 
F* = 70 and for densities (a) poa'^ = 0.07, (b) poa^ = 
0.20, and (c) poa^ — 0.70. The driving force is act- 
ing in y-direction as indicated by the arrow. The 
particle and box sizes in this figure are scaled to the 
same lengths for a better comparison of the structure. 



the range of the interaction, and k the reduced inverse 
screening length. In this work the energy is chosen to 
be Vb = lOfc^T, where fc^T is the thermal energy, and 
K — 4.0. The parameters are the same for the inter- 
action between particles of both types, and also for the 
cross interaction. The dynamics of the colloids is com- 
pletely overdamped Brownian motion. We neglect hydro- 
dynamic interactions which is justified for small colloidal 
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FIG. 2: The critical force Fc for lane formation versus 
density po- Brownian dynamics simulation results (cir- 
cles) are compared to the dynamic density functional the- 
ory (squares) . Solid and dashed lines are guide to the eye 
and suggest the phase boundary between states of 'no 
lanes' and 'lane formation'. The asterisks mark states of 
chosen densities (a) pocr^ = 0.07, (b) pocr^ — 0.20, and 
(c) poa'^ = 0.70 and a fixed force F* = 70 (indicated 
by the arrow), for which snapshots are shown in Fig. 1. 



volume fractions and long-ranged colloidal pair interac- 
tions. The constant external field Fi {i — +,~) is acting 
in y-direction of the simulation cell, shown in Fig.l, and 
drives particles of different types in opposite directions: 
F+ = F > for (-I-) particles and F^ = -F for (-) 
particles. We denote the dimensionless drive parameter 
by F* = Fa/ksT. In order to detect lane formation we 
apply the order parameter which is defined in Q . It is 
zero, when the system is isotropic, and one, when the col- 
loids are perfectly ordered in lanes throughout the whole 
system, each lane containing only particles of the same 
type. We could show in 4], that 0, averaged over many 
steady-state configurations, increases rapidly to values 
of about one when critical threshold value of the applied 
field, Fc is exceeded, indicating a nonequilibrium phase 
transition to the laning state. We define Fc to be the 
external field for which (/) — 0.5, implying that in av- 
erage half of the particles in the system are ordered in 
well defined lanes. Results for the critical force Fc ver- 
sus the total density are plotted in Fig. 2, showing a 
nonmonotonic behavior: for fixed high driving force (for 
instance F* = 70 as indicated by an arrow in Fig. 2) 
and increasing particle densities, there is first a transi- 
tion towards lane formation which is followed by another 
transition back to a state with no lanes. Corresponding 
snapshots for the three states at three different densi- 
ties poa"^ = 0.07,0.20,0.70 for a fixed force F* = 70 are 
shown in Fig. l(a)-(c). The average order parameter for 
these parameters is (j) = 0.11,0.80,0.05, respectively. At 
the highest density (Fig. 1(c)), there a still remnant of 
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laning visible but there are also regions of mutual jam- 
ming which destroy perfect lane formation. Hence there 
are qualitative differences in the no-lane phase at small 
and high densities. 

A dynamical density functional theory as described in 
Ref. accounts for the stabilization of the inhomogc- 
neous steady states beyond the critical drive. The density 
fields p± (r , t) obey the equation of motion which has the 
form of continuity equation: 



dp±{r,t) 

dt 



V • (f, t) + j'^> {r, t) + jt' (r , t)]. (2) 



T(2) 



•(3), 



The first current term j^^-* (r, t) is a diffusive current, and 
modeled as 0, • 



f±\r,t) 



P±{r,t)W f„ /-..N , (3) 



6p±{f,t) 



D = UbT/^ being the phenomenological diffusion coeffi- 
cient. J^[p+,p-\, the free energy cost for creating weak 
inhomogeneities around a homogeneous state, is given by 

= E j d\pu{r^\og[2p^{r)/p,] 
fc=± 

- \j (frd^r'c{\f - r'\)Ap{f)Ap{r') (4) 

with Ap(f) = p^(r) + /9_ (r) — po and c(r) the fluid direct 
correlation function, determined by V{r) [3. j^\f,t), 
is induced by the driving field in y direction. In a com- 
pletely demixed state, the field induces a Brownian drift 
velocity of Vd = F/j in the y direction, while in a com- 
pletely mixed state the flow is reduced due to mutual 
collisions between oppositely driven particle. These ef- 
fects can be incorporated via: 



Te„ — [1 
Po 



- ap±{f,t){pQ/2 



- p±{r,t))-l3p+{r,t)p^{r,t)], 



(5) 



where a is a phenomenological parameter to take care of 
the appropriate physical dimensions of the current and 
P takes care of the mutual collision effects. t) 
is also induced by the external drive but in the x di- 
rection transverse to the drive. Brownian collisions be- 
tween -|- and — particles having repulsive interaction will 
cause a displacement perpendicular to the field. The 
leading term contributing to j j_ (f , t) is the gradient in 

the density difference field l^j. j^\f,t) should scale 
with the total collision cross-section experienced per unit 
time. Let (Tq denote a typical range of the interaction 

— 1/2 

V{r) needed to perform a collision and as — Pq the 
mean interparticle separation. The number of collision 
event will be given by the number of particles of the two 
species in an area of a^crg- The typical time of colli- 
sion is ag/vd, so that the number of collision per unit 



time is ^p+{f,t){0)p-{f,t){as(To)'^. The cross-section of 
an individual collision is given by 7r(To/2, where the fac- 
tor of half takes care of the fact that the scattering is 
anisotropic due to the drive, only half portion of a parti- 
cle being exposed to collision. Consequently we find 

Ji''(f,i) = ±e,Tral/2—p+{f,t){0)p^(f,t){as(7of 

x^[p+(r,i)-p-(r-,t)]. (6) 

All currents remain unchanged under the interchange of 
the species, namely, p+(r, i) <-> p-{f,t) and F — > — _F, 
which satisfies the symmetry requirement, ctq has been 
estimated by an effective hard core diameter from the 
Barker-Henderson perturbation theory, (Tq = /□ [1 — 
ex-p{-V{r)/kBT)]dr. 

The dispersion relations for frequency uj with wave vec- 
tor q = {qxjQy)], obtained from the linearized equations 
of motion'l4l show that real frequencies, as required for 
the steady state bifurcations, can result for qy = 0, in- 
dependent of a and (3. The unstable wave vector in the 
y direction being qy = indicate structural homogeneity 
in the y-direction parallel to the drive as observed in the 
simulations. Two possible dispersion relations for qy ~ 
are: 



and 



where 



2F*ql-[l-po£{q.,Qv^0M (7) 



= -[1 - Pociqx-.qy = 0)]ql, (8) 



■"1,2 



^1,2^0 1 D denotes the dimensionless fre- 
quency and c{q) is the Fourier transform of c{r). 
remains negative for all qx, since 1 — pQc{q), being the 
inverse of the static structure factor, is a positive defi- 
nite quantity . w* changes sign, indicating the steady 
state bifurcation of a homogeneous state to an inhomo- 
geneous state [T^ . go where uj^ has a positive maximum 
as a function of g^;, will dominate the growth of the in- 
homogeneous phase. We approximate c(r) by that of an 
effective hard disk fluid with effective diameter (Tq for 
which analytical expressions are known We calcu- 

late numerically u}\ and locate graphically the instabil- 
ity point as well as go 01- ^^O'" ^ given po ^1 remains 
negative for all qx indicating the stability of the homo- 
geneous phase for very low F* . The homogeneous phase 
becomes unstable for F* > Fc to density perturbations 
with qx over a band about the maximum at qx — go- 
Note that a state with density modes having qy = but 
a finite qx = go indicates a lane phase as observed in 
the steady states of the BD simulations. The stability 
of the different phases has been shown over the F*-vs-po 
plane in Fig. 2, F* being marked by the squares for cho- 
sen densities.. The most striking feature of the stability 
diagram is the reentrant homogeneous phase as one in- 
creases the density for a fixed external force. Physically, 
the diffusive thermodynamic current tries to avoid any 
+/— interface via entropy of mixing, while the current 
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jj. (r, t) amplifies particle segregation in x direction via 
collisions induced by the drive. The diffusion becomes 
slower with increasing po that reduces the tendency of 
mixing. This requires stronger drive to ensure the colli- 
sion between two different species to have an appreciable 
jj. (r, t) contribution. This is what seen in the large po 
regime where Fc increases with pQ. However, the collision 
frequency itself goes down in the low po regime, and one 
needs stronger drive to counterbalance the fast diffusing 
current, leading to the reentrant behavior seen in Fig. 2. 

In conclusion, we predict a reentrance effect in binary 
mixtures of oppositely driven colloids for increasing den- 
sity. At fixed driving strength, there is a state with no 
lanes at small densities. This transforms into a state 
with lanes for intermediate densities and then there is 
reentrance of the state with no lanes. The reentrance 
behavior has a simple intuitive interpretation and should 
therefore be robust with respect to a variation of details 
of the interaction. We expect that it will show up also 
for asymmetric mixtures, in three spatial dimensions, and 
for different particle dynamics than the simple Brownian 
dynamics assumed in our model. In fact, recent simula- 
tions which include hydrodynamic interactions have 
revealed that the laning transition itself is indeed stable 



when explicit hydrodynamic interactions are added. 

The reentrance behavior should be verifiable in exper- 
iments. We think that the most promising candidate 
for an experimental test are binary mixtures of highly 
charged colloidal suspensions in a gravitational field. Re- 
cent progress has been made which highly deionized sus- 
pensions which can be subject to sedimentation at 
different colloidal densities. Another possibility are mix- 
tures of colloids and polymers under gravity |2ll |. The 
reentrance effect may finally provide also a quantitative 
framework to explain different states of pedestrian dy- 
namics in pedestrian zones for varied population density 
[2^ : The low-density state without lanes is typical for 
low populations. For intermediate densities, pedestrians 
form lanes, but for even higher densities the lanes are 
broken again and there is a crossover to panic dynamics. 
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